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We analyze various scenarios for entangling two initially unentangled qubits. In particular, we 
propose an optimal universal entangler which entangles a qubit in unknown state with a qubit 
in a reference (known) state |0). That is, our entangler generates the output state which is as close 
as possible to the pure (symmetrized) state (|^I')|0) + |0)| >]/)). The most attractive feature of this 
entangling machine, is that the fidelity of its performance (i.e. the distance between the output and 
the ideally entangled - symmetrized state) does not depend on the input and takes the constant 
value T = (9 + 3v2)/14 — 0.946. We also analyze how to optimally generate from a single qubit 
initially prepared in an unknown state |$) a two qubit entangled system which is as close as possible 
to a Bell state (|*)|* ± ) + |* ± )|*)), where (*|* ± ) = 0. 
PACS number: 03.67.-a, 03.65.Bz 
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I. INTRODUCTION 

A pure quantum state of two systems A and B is said 
to be entangled if it is not a product of a state for A 
and a state for B. Two systems in an entangled state are 
correlated, and these correlations are intrinsically quan- 
tum mechanical M. For example, one must use entangled 
states in order to produce violations of Bell inequalities 
or in the test of local realism proposed by Hardy [§,1). 
Entangled states also play a key role in quantum infor- 
mation, in particular they are essential in quantum tele- 
portation Q and in superdense coding In quantum 
computers entanglement is one of the features of quan- 
tum mechanics which give these machines their power 
§• 

Here we would like to consider the problem of how to 
produce entanglement. In particular, if we are given par- 
ticles, or systems, A and B in the pure states \^) A and 
| $) B , respectively, we would like to produce the state 
(|*) A |3>) B + |$)aI*)b) (up to normalization). Formally, 
we are looking for the symmetrization map 



(1.1) 



In what follows, where possible we omit explicit sub- 
scripts A and B. The order in which the vectors are 
written in the tensor products implicitly denotes to which 
system they belong (i.e. the left vector corresponds to the 
system A, while the right vector corresponds to the sys- 
tem B). We assume that the two quantum systems (e.g., 
qubits) are physically distinguishable. For instance they 
could be located in different regions of space. The task 
is to entangle their internal degrees of freedom. 

That the symmetrization cannot be done perfectly via 
a unitary transformation can be shown by the following 
argument. We consider the case in which \9) and |$) 
are both qubits. A perfect transformation would have to 
transform the basis vectors as 



|oo)K) 
|oi)N) 

|io>N) 
|n)N) 



t'i) 



_(|01> + |10»h) 

-L(|oi) + |io))|« 3 ) 
|ii>K>, 



(1.2) 



where the \vj), for j = 0,4, are normalized "machine" 
vectors, i. e. we assume that the entangler itself has its 
own degrees of freedom. In addition, it is assumed that 
the entangler is always initially in the same state, \vq). 
Unitarity requires that (W2IW3) = 0. Now let us consider 
the case where the input vectors are \^f) = a|0)+/3|l) and 
1$) = |0) (i.e. the state of the qubit A is unknown, while 
the qubit B is in a known state). The transformation 
(|1.2[) gives us 



|tf)|0) ->a|00>|«i) 



V2 



(|01) + |10))h), (1.3) 



whereas what it should produce is a vector proportional 
to |#)|0) + |0)|#), which in the basis |0), |1) reads 

|*>|0) |*)|0> + |0>|¥>=2a|00>+/9(|01> + |10>). 

(1.4) 



The vectors in the right-hand sides of Eqs.(1.3) and (1.4) 
are clearly not the same, no matter what choice is made 
for and |i>3). Therefore, we need to search for devices 
which will produce approximate versions of the desired 
state or will produce this state but with a probability 
which is less than one. 

One way of creating a symmetrized state out of two in- 
dependent systems is by means of a measurement - that is 
the two systems are optimally measured and their states 
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are estimated. Based on this estimation a two-particle 
entangled state is prepared. If we begin with two qubits 
prepared so that one of the states is known (|0)) while the 
other is unknown (| I J / )), we need only estimate the state 
of one of the particles and this can be performed with a 
fidelity equal to 2/3 The information gained from 

the optimal measurement is then used i n the preparation 
procedure. This is discussed in Section II A. 

We shall present quantum mechanical entangling 
transformations which generate entangled states with 
much higher fidelity than can be achieved by measur- 
ing the input particles. In Section [I B we briefly discuss 



a probabilistic symmetrization (entanglement) which can 
be realized via a controlled-SWAP gate. The probability 
of success in this procedure is input-state dependent. In 
Section [II we present the optimal input-state indepen- 
dent quantum entangler and we also study the insepara- 



bity of the outputs of this entangler. In Section |V| we 
show that the universal-NOT gate |l can also serve as a 
very interesting entangling device. 



II. STATE-DEPENDENT SYMMETRIZATION 

We shall first look at two examples of processes which 
produce entangled states, for which the quality of the out- 
put depends on the input state. That is, these procedures 
work better for some states than for others. The first is 
perhaps the most obvious method, we simply measure 
the input state. We shall consider a more limited prob- 
lem in this case, entangling an unknown with a known 
state. The output state resulting from this procedure is 
only an approximation to the desired one. The second is 
a probabilistic method; the output when it is produced 
is ideal, but the probability of successfully producing it 
is less than one. In this case we shall consider the full 
problem of entangling two unknown states. 



output will be determined by calculating the fidelity be- 
tween the actual output and the desired output. We shall 
first examine a specific strategy and then find an upper 
bound on the fidelity for a wide class of measurement- 
based procedures. 

Our first measurement-based scenario can then be re- 
alized in the following way. In the case of a single in- 
put qubit the optimal way to estimate the state, is to 
measure it along a randomly chosen direction in the 
two-dimensional Hilbert space @J|. Therefore, the first 
step in implementing the measurement-based procedure 
is choosing a random vector (77), where 



|r ? )=cosy|0)+e^' sin^|l), 



(2.2) 



and measuring \*5f) along it. If the result is positive, then 
the output is taken to be |$) AB , and if negative, the out- 
put is \&} AB , where 



\ o)\v) + \v)\Q) 

2(l + cos 2 ^ 



2cos^|00) + \/2V v ' sin f |+) 
2(1 + cos 2 f) 



(2.3) 



and 



|^}|0) 



2(1 



sin 2 f) 



2e- i * / sin 00) 



- V2cos^ 



2 1 



sin 



2 i3' 



(2.4) 



where the state {r/ 1 -) is the state orthogonal to \rj), 







d' 



= QT lLp sin— 10) - cos— 11 



(2.5) 



A. Entanglement via measurement 

Our task is to entangle an input qubit in an un- 
known state with a reference qubit in a known state 
|0). That is, we want to realize the symmetrization map 
|0) A |4 , ) B — > l^w) AB with the output parameterized as 



|*(«) )ab = 



2cos||00) 



-V2e^ sinf |+) 



/2(l + cos 2 f) 



(2.1) 



The approach we will discuss here is as follows: firstly, the 
unknown single-qubit state | "J) is measured and then us- 
ing the information gained thereby an approximate ver- 
sion of the desired output is constructed. In order to 
specify this procedure in more detail, we must describe 
what measurement is to be made and how its results will 
be used to construct the output state. The quality of the 



For a particular orientation of the measurement appa- 
ratus, i.e. for the particular choice of the state (77) this 
measurement-based scenario gives the two-qubit output 
density matrix 

^«*)(^|<?V) = K^ir?)! 2 !*)^! + 1(*|^)| 2 |$)($|. 

(2.6) 

To get the final output density matrix one averages this 
over all possible choices of the measurement (i.e. over all 
vectors |?7)) 



47T 



(2.7) 



Finally, the fidelity can be found by computing the ma- 
trix element of this density matrix in the ideal output 
state, |*( id )).„, 



2 



F{Q,<p) = (^ d) \p {out) (ti,<p)\y (ld) ) 



(2- 



This fidelity depends on the input state, and this de- 
pendence can be eliminated if we average over all input 
states 



T 



(2.9) 



This is the proper fidelity to use to judge how well our 
proposed strategy performs if we assume that all input 
states are equally probable. A more explicit expression 
for it is 



T 



1 



271 



2 71 



dtp / dip' 



sin ■Odd 



167T j j j Jo 

\( v \*)\ 2 IW)| 2 + |(^|*)| 2 1<*|*)| 2 ' 



Explicitly evaluating this integral we find 
T = 54 + 112 (In 2) 2 - 154.5 In 2 ~ 0.719 



sini?W 

(2.10) 

(2.11) 



which is a bit larger than 2/3, the fidelity of the estima- 
tion of a state of a single qubit. 

Let us now generalize this procedure. We shall again 
begin by choosing a random vector but now according 
to a distribution which we shall leave unspeci- 

fied for now. The output density matrix is taken to be 
either p\ (rj) if the measurement result is positive or po (77) 
if it is negative, where 

pM = J dn" Pj (r,^\di,p')\T(r 7 p"))(r(r,p")\, 

(2.12) 

with j = 0, 1, and 

|r(tf'V)) AB =co S y|00)^ B + e ^"sin^|+)^ B . (2.13) 

The conditional probabilities pj will also be left un- 
specified; this allows us to consider a wide class of 
measurement-based strategies. The output density ma- 
trix, for a particular I77) is then 



Averaging over \r)} gives us the final output density ma- 
trix 



(2.15) 



and the fidelities for a specific input stat e an d aver age d 
over all input states are given by Eqs. (2.8) and (|2.9]) , 
respectively, bu t w ith py° ut ) computed from Eq. ( |2.15| ) 
instead of Eq. (2.7). In particular we have that 



t= / / dtfdrt'Y,fj&'yrfW)Pj(f',</';*',<//), (2-16) 

J J j=0 
where 

p j (r,p"-,'d',p')=p j (r,p'Y,p') q ^',p'), (2.i7) 

for j = 0, 1, and 

f 1 .0 1?" 

/o = j / ^ 2(l + cos 2W 2)) |2COS 2 C ° S T 
+V2e 1 ^"-^ sin - sin — | 2 | (*|^) | 2 



A = / dn 



1 . 1? 1?" 

-2(1WW2» |2COS 2 COS T 
+V2e l(v "-^ sin t gin 1L | 2 | (V\r)) | 2 . 



(2.18) 



p(v) = \(vm 2 pi(v) + \(v ± m 2 po(v)- 



(2.14) 



What we can now do is to find an upper bound for the 
fidelity, J 7 , for any distribution of the vector \rj) and any 
prescription for using the result of the measurement along 
\rj) to manufacture the entangled state. We note that for 
3=0,1 

1 = J d£l' fdn"P j ('&", </?";i?', <ff), (2.19) 

which implies that 

F^supl/ol+suplAI, (2.20) 

where the supremums are taken over the range < 
1?',1?" < 7T and < <p',<p" < 2tt. 

Our first task is to find explicit expressions for the 
functions fo and f\. We have that 



fo = di cos — sin — 



, 2 2 & 

a? cos — cos — 
2 2 



1 , 2 1?" 2 1}' 
— as sin — cos — 

2 2 2 



/i = 



2 i?" 2 i?' , „0" J' 
d\ cos — cos — + d-2 cos — sin — 



1 j • 2 0" • 2 1?' 
—do sin — sin — 

2 2 2 

/s, / // a ^" ^ • 0" • 0" 
V 2(22 COS(<£ — ^3 ) cos — cos — sin — sm — 

1?' 



1 , . 2^" 

-efesm T cos - 



1 , • 2 ^" ■ 2 & fK, / // A tf" ^ ■ • 

--fl3sm — sm h v 2a2 cos(ip — 95 ) cos — cos — sm — sm — 

2 2 2 2 2 2 2 



(2.21) 
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where 



dx = 2 In 2 - 1 
d 2 = 3 - 4 In 2 
d 3 = 81n2-5. 



(2.22) 



From the above equations it is clear that in order to max- 
imize /o we need to choose ip" — tp' = n and to maximize 
/i we need to choose ip" — ip' = 0. Making these choices 
and simplifying the resulting expressions we find that 

/o(tf",7r;tf',0) = -[1 + d costf" - c 2 cos 0" cos i?' 

+C3 sin 1?" sin 1? ] 
/i(i?",0;i?',0) = -[l + c x cosi9" + c 2 costf"cosz9' 



where 



+c 3 sini?" sin &'], 



Cl = 3 - 41n2 

c 2 = 121n2 - 8 

c 3 = V2(3-41n2). 



(2.23) 



(2.24) 



These functions can now be maximized. The maximum 
of /o occurs at 1?' = tt and <&" = 0, and the maximum 
of fi occurs when d 1 — and 1?" = 0. The maximum 
values of both functions are the same and are approxi- 
mately equal to 0.386. This implies that the fidelity for 
this kind of a measurement-based strategy must satisfy 



T < 41n2- 2 ^ 0.773. 



(2.25) 



As we shall see, a method which maintains quantum co- 
herences at all stages of the process can do better than 
this. 



B. Controlled-SWAP gate 

We now begin with systems A and B of the same phys- 
ical origin. Their pure states are described by vectors in 
the D-dimensional Hilbert space Tt, so that both together 
are described by Tt <g> Tt. Let = 1, ■ ■ ■ D} be an 

orthonormal basis for Tt. System A is in the state 



D 
3 = 1 

and system B is in the state 

D 



(2.26) 



l*) B =E d 



(2.27) 



Our objective i s to produce the (entangled) symmetrized 
state [see Eq. (1.1)] 



D D 

+ mm =^^4+^)1^)1^). (2.28) 
j=i k=i 

(here we omit the normalization factor). 

Recently Barenco et al. [[[o| have shown that the en- 
tanglement (symmetrization) of the form ([0]) can be 
performed when the two input qubits interact via a 
controlled- SWAP (Fredkin) gate with an ancilla initially 
prepared in a specific state. The entanglement is achieved 
when a conditional measurement is performed on the an- 
cilla. Exactly the same scenario can be used not only for 
qubits but for arbitrary quantum systems. To show this 
we briefly review the operation of the controlled-SWAP 
gate. 

This gate has three inputs. The first, the control bit, is 
a qubit. The second and third are for Z?-dimensional sys- 
tems. The control bit is unaffected by the action of the 
gate. If the control bit is |0), then the gate does nothing, 
i.e. the output state is the same as the input state. If 
the control bit is |1), then the two D-dimensional states 
are swapped. This can be accomplished by the following 
explicit unitary transformation: 

|o)K)K) - |o)K)K> ; 
|i)K-)K)-|i)K)K)- (2.29) 

Summarizing, the action of our controlled-SWAP gate is, 



|0>|¥>|$>->|0>|*>|<l»> 
|1)|¥>|$>->|1>|<I»>|¥) 

We now define the qubit states 



(2.30) 



|« + > = i(|0) + |l»; |«_) = -L(|0> - |1», (2.31) 

and take the input state of the controlled-SWAP g ate to 
be \v+)m A m B . Using the SWAP transformation Q2.3C| ) 
we find that the output state is 

l^) ) = _L ( | 0) |^ ) | $) + | 1) | $) ^ )) 

= i|« + >(|*>|*> + |3>|*>) (2.32) 

+ i|«_)(l*>l*>-|*>l*». 

If we now measure the qubit in the \v±) basis we ob- 
tain the states (|*)|$) ± I*)!*)) with probabilities (1 ± 
|(\j/|$)| 2 )/2, respectively. As we see the probability of 
generation of a particular entangled state explicitly de- 
pends on the (unknown) states of the two systems. In 
particular, let us assume we begin with two orthogonal 
qubits, |\&) and l^ -1 ). Then either of the maximally en- 
tangled state,(|*)|*- L )±|*- L )|'I'))/v / 2 can prepared with 
probability 1/2. 

We stress that the probability of the success in this 
entanglement (symmetrization) procedure is input-state 
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dependent. In what follows our task will be to find a 
"machine" which entangles the input with a constant (i.e. 
input-state independent) fidelity. This covariance prop- 
erty of the entangler with respect to unitary transforma- 
tions performed on the input qubits makes the entangler 
universal. 



III. UNIVERSAL ENTANGLERS 



where |±) are symmetric and anti-symmetric Bell states 
in the given basis 

|±> = -J=(|01>±|10». (3.6) 

In what follows we will briefly discuss the anti-symmetric 
entangler and then we will concentrate on the symmetric 
entangler. 



Suppose we again consider the problem of constructing 
a device which will entangle a qubit in an arbitrary un- 
known state \^f) = a\0) + (3\1) with a qubit in a known, 
reference state, which we shall take to be the basis state 
|0). Before we proceed further we have to specify proper- 
ties of the entangling map. In fact, we can consider two 
maps. The symmetrization map 



S:\0) A m 



l* H) L B =^(W|o) + |o)|*)), (3.i) 



and the anti-symmetrization map 

= JV a (|¥)|0)-|0)|¥)), (3.2) 



A:\0) A \*) B -+\^ d) ) 



where N a , s are corresponding normalization factors. As 
we have shown in the introduction perfect entanglers for 
arbitrary unknown states cannot be constructed. So 
the task of the physically realizable symmetric (anti- 
symmetric) entangler is to produce outputs as close 
as possible to the ideally entangled states |^^^)^ B 

(I*' ') AB ). In what follows we will quantify the qual- 
ity of the performance of the universal entangler with 
the help of the fidelity 



T 



(\fr( id ) \p( out ) \\fr( id )j 



(3.3) 



We shall impose the condition that the value of th is fi - 
delity does not depend on the input. The fidelity ( |3.3j ) 
is a good measure of the accuracy with which the en- 
tangler produces the desired output state, but we would 
also like to evaluate the degree of entanglement of the 
actual output state. Here, however, we have a problem 
which is due to the fact that it is still not clear how to 
quantify the entanglement of a quantum system which is 
in a mixed state. When a bipartite system is in a pure 
state, then the von Neumann entropy of subsystems can 
serve as a measure of entanglement. In the case of impure 
states more sophisticated measures are required (see for 
instance [0-^3)). 

In terms of the basis vectors, the input state is a|00) + 
/3|01), and the ideal output state in the case of sym- 
metrization is 



(2a|00) + V2/3|+)) 
1 ' (4|a| 2 + 2|/3| 2 ) 1 /2 ' 



(3.4) 



while in the case of the anti-symmetrization we have 

= |-), (3.5) 



A. Entanglement via anti-symmetrization 

Recently Alber @ studied a quantum entangler which 
takes as an input a quantum-mechanical system prepared 
in an unknown pure state \^) A and a reference (known) 
state (let us say \0) A ) and at the output generates a two 
particle entangled state p A g^ which is optimally entan- 
gled. Alber imposed two constraints on the output of the 
universal quantum entangler 



Tr 



(out) 



= Tr t 



(out) 



1 

D 



and 



(out) 

" AB 



minimum. 



(3.7) 



(3.8) 



Where D is the dimensionality of the Hilbert space of 
the system A (B) and S is the von Neumann entropy 
S = — Trpln p associated with a given density operator p. 
The first condition corresponds to the requirement that 
the subsystems at the output are in the maximally mixed 
state while the second conditions guarantees that the 
whole system is as close as possible to a pure two-particle 
state. Alber has found the solution for this problem. It 
turns out that the two-particle state which is produced by 
the optimal (with respect to the above conditions), uni- 
versal entangler is independent of the input state and 
is equal to a maximally disordered mixture of all possible 
anti-symmetric Bell states. In the case of qubits (D = 2) 
there is only one possible anti-symmetric Bell state |— ). 
That is, Alber's machine realizes the anti-symmetric en- 
tangler. We see that the universality of Alber's entangler 
means that all inputs are mapped to a single output (the 
anti-symmetric Bell state |— )), so the ideal output state 
is a priori known, and one could instead build a device 
which just prepares the known output state. In the anti- 
symmetric entangler the information initially encoded in 
the qubit A is completely lost. But our task is different, 
we want to redistribute the initial unknown information 
encoded in the state of the qubit A, into the entangled 
state of two qubits. Therefore we will analyze universal 
entanglement via symmetrization, because the ideal state 
( |3.4[ ) directly contains information about the initial state 
of the qubit A. In other words, we consider the entan- 
gling procedure not only as the way to generate the state 
with highest possible entanglement but also we require 
that this state contains as much information about the 
input(s) as possible. 
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B. Entanglement via symmetrization 

Let us now construct a machine which entangles an 
unknown state with the known state |0). Taking into ac- 
count the basic features of the symmetrization transfor- 
mation (3.1) we can assume that the basis vectors trans- 
form as 



|oo>h> 
|oi>K) 



|oo)K) 
|oo>K) 



-)ko) 



(3.9) 



where |u>o), \wi), \xq), and |a;i) are states of the entan- 
gler itself. The entangler is initially always prepared in 
the state \v ). 

We want to impose the condition that the fidelity be- 
tween the actual output state and the ideal output state 
be independent of the state \^}, but before doing so let 
us state the restrictions which unitarity places on the 
machine vectors. These are 



INI 
(wo\wi) 



x \r = l; 
2 = i; 
) = 0, 



+■ IFill 
<ar |ari 



(3.10) 



where ||a;|| 2 = (x\x). We now calculate the output two- 
qubit de nsit y matrix p( out ) by using the transformation 
in Eq. ( |3.9| ) to find the full output density matrix and 
then tracing out the mac hine degrees of freedom. We 
then find the fidelity (3.3) by taking the matrix element 
of this density matrix in the ideal output state. Our task 
is to find the machine vectors \x -) and \wj) (j = 0, 1) 
such that the fidelity T does not depend on the input 
state |\&) and simultaneously is as close as possible to 
unity. 

We find that if we choose |xo) to be orthogonal to each 
of the other machine vectors and \w\) to be orthogonal 
to |xo) and |i«o), then the output fidelity will be inde- 
pendent of the phases of a and j3. Making these choices 
we find that 

F=N- 1 {2\a\*\\w f + \f3\*\\x x f + \a\ 2 \f3\ 2 

V2((w \ Xl ) + ( Xl \w )) + 2\\ Wl \\ 2 + ||x || 2 ] } , (3.11) 

where N = 2\a\ 2 + \fi\ 2 . 

In order for this expression to be independent of \a\ 
and \f3\ it is necessary that the expression in the curly 
brackets be proportional to 



(2M 2 + l/?l 2 )(k 



J ) = 2|a| 4 + 3|cf 



(3.12) 



Comparing this expression to Eq. ( 3.11 ) we see that 



as large as possible. If we now make use of the unitarity 
conditions and the two equations above, we find that 



\ w o\ 



2 r 1 
1 - -V2cos^ = — 

3 Ikoll 2 



where 



cos 



_ (xi\w ) + (w \xi) 



(3.14) 



(3.15) 



From Eq. (3.14) we see that ||wo|| 2 will be a maximum 
when cos p = 1, which implies that |iuo) and \xi) are 
parallel. When this condition is satisfied, we find that 



J 7 =\\w \\ 2 = 



9 + 3\/2 



14 



(3.16) 



which gives 0.946 as the approximate value of the fidelity. 
This means that the output state p( out > is indeed very 
close to the ideal state, and it should be remembered 
that this fidelity is the same for all input states. 

We can summarize our results for the machine vectors 
as follows. From the above analysis we see that we can 
take the machine state space to be three dimensional. 
Define 



cos( 



9 + 3V2 
14 



1/2 i- -i 1/2 



smt 



5-3\/2 



14 



(3.17) 

and let {\vj) \j = 1, . . . 3} be an orthonormal basis for the 
machine vector space. We then have 

|?«o) = cos 6\vx) 
\wi) = sm9\v 2 ) 

\xo) =sm6\v 3 ) (3.18) 
|aci) = cos^|ui), 
and our transformation in terms of basis vectors becomes 

|00) |u ) -» cos 0\OO)\vi) + sin 0\+) \v 3 ); 

|01)|«o) -^sin6»|00)|w 2 ) +cos0|+)|ui). (3.19) 

By construction this is the optimal entangling transfor- 
mation which entangles an unknown pure state with a 
known reference state. 

Alternatively, for = a|0) +/3|1) we can rewrite this 
transformation in the form 

|O)|*)|u o )^cos0(a|OO) + /3|+))K) 

+ sm0(a\+)\v 3 ) + /3\00)\v 2 )) . (3.20) 

When the trace over the entangler is performed we ob- 
tain the density operator p^°^ describing the two qubits 
A and B at the output of the quantum entangler 



\\ w o\\ = iFill 

3|| Wo || 2 = V2({xi\w ) + {wolxt}) + 2\\ Wl \\ 2 + \\x \\ 2 . (3.13) 

If these conditions are satisfied, then the fidelity is simply 
equal to || wo || 2 , so that we want to make this quantity 



\a\ 2 cos 2 ( 



\a\ 2 sin 2 i 



f sin 2 9) 1 00) (00 1 
| 2 cos 2 0)|+)(+| (3.21) 



(out} 

Pab 



It is important to stres s tha t the fidelity ( |3.3| ) associated 
with the output state (3.20) is input state independent. 



cos 2 (9(a/T|00)(+| + (00|) 



6 



C. Remarks 



0.7 



Throughout this paper we have utilized the fidelity 
( |3.3D as the measure of the performance of the quantum 
entangler. The universality (covariance) of the entangler 
is expressed in the fact that the value of the fidelity T 
is equal for all input states. We note that this covari- 
ance constraint is equivalent to the requirement that the 
Bures distance [fl5| defined as 



dn(pi,/fc) =\/2(l-Tr 



1/2. .1/2 
Pi P2Pl 



1/2 



(3.22) 



between the ideal state and the output of the en- 

tangler p A °^ is constant. In our particular case we find 
the Bures distance to be 



d B = 2sin(0/2) ~ 0.0541, 



(3.23) 



for all inputs. This distance is very small indeed. It is 
important to note that the Hilbert-Schmidt norm 



2-|!/2 



dHs{pi,P2) = [Tr(pi - p 2 
which in our case can be expressed as 



l-2^ + Tr(p( l o ; t ))' 



1/2 



(3.24) 



(3.25) 



0.6 
. 5 
0.4 
0.3 
0.2 
0.1 




a 2 



FIG. 1. The von Neumann entropy of the single-qubit 
state p A when the two-qubit system is in an ideally entangled 
state |\j/^2) (line 1) and when the output state p^ J) is given 



by Eq.( 3,21 ) (line 2). In both cases we assume a and j3 to be 
real. 

It is interesting to find the entropy of the two-particle 
state pVg*' at the output of the entangler as a function of 
the initial state (in the ideal case the two-particle system 
is always considered to be in a pure state with 5 = 0). 
We plot this entropy in Fig. |^. We see that the total 
entropy of the output is state-dependent and it takes the 
minimal value for a 2 = 1/2. Therefore the entropy of the 
subsystems does not indicate whether they are entangled. 



is not input-state independent because Tr(p^^) de- 
pends on the initial state. This is closely related to the 
fact that the von Neumann entropy of the state p A °^ is 
state dependent (see below). 



D. Inseparability of the output qubits 

We note that the entanglement between the two qubits 
prepared in the state depends on the particular 

form of the state |<F) = a|0) + /3|1). Because \^^) 
is a pure state we can quantify the degree of entan- 
glement via the von Neumann entropy S of one of the 
two qubits under consideration, i.e. Sa — — TrL^lnp^] 
(obviously Sa — Sb)- For a = 1 the entropy is equal 
to zero, which corresponds to a completely disentangled 
state (we note that in this case = |0)|0)). The 

entropy takes the maximal value S = In 2 for a = 
when = (|0)|1) + |l)|0))/%/2. We plot this en- 

tropy in Fig. [l] (see line 1). The entropy of the individ- 
ual particle (qubit) at the output of the entangler, i.e. 
p{out) _ rp r ^(oMt) j g a j wa y S l ar g e r than in the ideal case 

(see line 2 in Fig. |l]). Nevertheless, for the case a = 
we have in this case S(a = Q) = 0.998 In 2, i.e. this en- 
tropy is very close to the entropy of a qubit in the ideal 
case. Unfortunately, this entropy in the case of an im- 
pure two-particle state cannot be used as a measure of 
entanglement. 




a 



FIG. 2. The von Neumann entropy of the two- qubit state 
Pab^ a t t ne output of the entangler [see Eq. (3.21)] as a func- 
tion of a 2 . We assume a and /3 to be real. 

We need to check whether the two qubits A and B 
at the output are indeed quantum-mechanically entan- 
gled. Quantum-mechanical entanglement of two qubits 
formally means that the density operator of these two 
qubits is represented by an inseparable matrix (s ee M ) 
It follows from the Peres-Horodecki theorem that |l6p7| 
the necessary and sufficient condition of inseparability 
of the two-qubit density matrix p AB is that the corre- 
sponding partially transposed matrix p^ 2 has at least 
one negative eigenvalue. 
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For instance, let us consider the state |^) = a\0) + 
with real amplitudes a and j3. The partially trans- 
pos ed m atrix corresponding to the state given by 
Eq.(3.4) has one negative eigenvalue 



E(a) 



2(a 2 



1 



(3.26) 



We plot this eigenvalue in Fig. || (see line 1). We see 
that the eigenvalue is negative for all values of a except 
a = 1 when | = |0)|0). The minimal value of the 

eigenvalue is achieved for a = when the two qubits are 
in the maximally entangled state ( 1 01) + \10})/y/2. 



-0 . 3 




a" 

FIG. 3. Here we plot the negative eigenvalue Eq. j3.26| ) of 
the partially transposed matrix of the density operator p { ^ eal ^ 
when the state |^} has real amplitudes a and /3 (see line 1). 
The negative eigenvalue of the partially transposed matr ix as - 



sociated with the density operator p^g given by Eq. (3.21) 



as functions of a is presented by line 2. (We assume a and P 
to be real.) 

Now we utilize the Peres- Horodecki theo rem to check 
whether the state pv£*' given by Eq. ( 3.21 ) describes an 
entangled state of two qubits. Firstly, we find that the 
partially transposed matrix corresponding to the density 
operator (3.21) has one eigenvalue which is negative for 
all values of a (here we assume a and (3 to be real). In 
particular, this eigenvalue for a = is 



E(a 



cos v — 



e - (cos 4 e + sin 4 e) 



1/2 



(3.27) 



which is the minimal value (~ —0.447) of the nega- 
tive eigenvalue. On the other hand the maximal value 
(~ —0.001) is attained for a = 1 



E(a = 1) 



sin 2 (9- (cos 4 + sin 4 6 



1/2' 



(3.28) 



The complete dependence of E(a) is shown in Fig. [|. 
From this figure we clearly see that the output density 
operator is inseparable for an arbitrary input considered 
in this Section. We note, that if the entanglement is 
measured in terms of the tangle as introduced by Woot- 
tcrs |13] then the negative eigenvalues E of the partially 



transposed density operators perfectly reflect the degree 
of entanglement between the two qubits in our cases. 

By construction the fidelity of the entangler in this 
case is constant but the actual degree of entanglement 
is state-dependent. This suggests that it would be inter- 
esting to find an entangler, whose output states have the 
same degree of entanglement irrespective of the input, 
yet still carry information about the input. 



IV. ENTANGLEMENT VIA UNIVERSAL NOT 
GATE 

Even though the negative eigenvalue of the partially 
transposed density matrix cannot be directly used as the 
measure of entanglement, we see that the degree of en- 
tanglement between two qubits generated in the entan- 
gler (3.19) depends on the input state. In what follows 
we describe a different type of the entangler, which out 
of a single qubit \^f) generates a two-qubit state as close 
as possible to the state 



I*) 



|{*,* X }> = (|*)|^) + \^)\^))/V2. (4.1) 



We will present an entangler which not only produces 
the state which is as close as possible to the ideal state 
}) but also has the property that the fidelity does 
not depend on the input state. In addition, the degree 
of entanglement also does not depend on the input. This 
type of the entangler implicitly assumes creation of the 
state l^" 1 ") from the input That is, we face the prob- 
lem of creating an orthogonal state from unknown input. 

It is not a problem to complement a classical bit, i.e. to 
change the value of a bit, a to a 1 and vice versa. This 
is accomplished by a NOT gate. Complementing a qubit, 
however, is another matter. The complement of a qubit 
l^) is the qubit l^" 1 ) which is orthogonal to it. But it is 
not possible to build a device which will take an arbitrary 
(unknown) qubit and transform it into the qubit orthog- 
onal to it. As shown in Ref. || the ideal universal-NOT 
(U-NOT) operation corresponds to the inversion of the 
Block (Poincare) sphere. This inversion preserves angles 
(related in a simple way to the scalar product of 
rays), so by Wigner's Theorem the ideal U-NOT must 
be implemented either by a unitary or by an anti-unitary 
operation. Unitary operations correspond to proper ro- 
tations of the Poincare sphere, whereas anti-unitary op- 
erations correspond to orthogonal transformations with 
determinant — 1 . Clearly, the U-NOT operation is of the 
latter kind, and an anti-unitary operator O (unique up 
to a phase) implementing it is 



e(a|0> =/3*|0> 



a* 1) 



(4.2) 



The difficulty with anti-unitarily implemented symme- 
tries is that they are not completely positive, i.e., they 
cannot be applied to a small system, leaving the rest of 
the world alone. 
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Because we cannot design a perfect Universal-NOT 
gate, we have introduced in Ref. j9j an approximate opti- 
mal U-NOT gate (an analogous spin- flip operation has re- 
cently been introduced by Gisin and Popescu Eq] ) . This 
device takes as an input the qubit A in the state and 
generates at the output a qubit in a mixed state as close 
as possible to the orthogonal state I 1 ?" 1 "). The role of the 
U-NOT gate is played by two additional (ancilla) qubits 
B and C. So, all together the transformation involves 
three qubits and it can be explicitly written as 



\*)AX) 



AS I IC 



+ 71 |{#,¥- L }) 



I*) 



(4.3) 



where \X) BC , is the initial state of the U-NOT gate; 
70 = \/2/3 and 71 = — \/l/3. In this particular transfor- 
mation the qubit C at the output is in the state which is 
as orthogonal as possible to the input state. The fidelity 
of this transformation is input-state independent and is 
equal to T = 2/3. 



A. U-NOT as the entangler 

It is interesting to note that the two -qubit state p^°^ 
at the output of the U-NOT gate ( |4.3| ) has the form 

P [ :f =7il{*,* ± })({*,* ± }|+7ol^)(^l- (4-4) 

The me an fi delity between the state pv^' and the ideal 
output J4.l| ) is input-state independent and takes the 
value T = 1/3. This again corresponds to the fact that 
the Bures distance between the actual output of the en- 
tangler and the ideal output is input state independent 
and equal to oIb = (2 — 2/v / 3) 1 ^ 2 - We can easily check 
that the partially tra nspo sed matrix corresponding to 
the density operator ( |4.4| ) has one negative eigenvalue 
E = (2 — v5) /6 which is constant and does not depend 
on the initial input state 

We note that the Universal NOT gate ( [4.3[ ) acts also a 
quantum doner, i.e. the two qubits A and B are the op- 
timal clones of the input (for details see Refs. ||@). It 
is the optimality of the transformation (fll|) with respect 
to cloning and the generation of the optimally orthogo- 
nal state (i.e. the univer sal NOT gate) which indicates 
that the transformation (4.3) also serves as the optimal 
universal entangler. 



(perfect) entangler is due to the linearity of quantum me- 
chanics. On the other hand, another consequence of the 
linearity of quantum mechanics is the fact that the en- 
tangled quantum-mechanical states cannot be used for 
super-luminal communication. Gisin pl| has shown that 
this no-signaling constraint implies bounds on the fidelity 
of universal cloning and the universal U-NOT gate. In 
the case of cloning the bound on fidelity is T = 5/6, while 
in the case of the U-NOT gate the bound is T = 2/3. We 
note that the transformation (4.3) achieves both these 
bounds when used as the doner or the U-NOT gate, re- 
spectively. Recently Alber Jl4|] used this idea of Gisin to 
prove that the upper bound in the fidelity of the anti- 
symmetric entangling is equal to unity. The no-signaling 
constraint can also be used to derive an upper boun d on 
the fidelity of the entangling operation given in Eq. (4.1) 
p2f . We will present a proof, which is based on the meth- 
ods developed in reference pjj , that this upp er b ound is 
T = 1/3, which means that the U-NOT gate (4.3) serves 
as t he optimal universal entangler in the sense of Eq. 

We consider a process in which a single particle input 
state is mapped into a two particle output state. The 
input state can be represented as 



P {m) {rh) = ^(1 + fh-B) 



1 



(4.5) 



where to is a real vector whose length is less than or 
equal to unity. The most general two-particle output 
state, which is hermitian and has a trace equal to one, 
can be expressed as 



(out) 



1 



(to) = -\t + a- B <g> 1 + 1 <g) 6 • B 

j,k=x,y,z 



(4.6) 



where a, b, and tjj~ are functions of to. The requirement 
that the reduced density matrixes of the two output par- 
ticles be the same, which we shall impose, implies that 
a = b. 

We now want to impose the requirement of covari- 
ance. This means that if p( m \m) is mapped onto 
p\° ut )(m), and if u is a matrix in SU(2), then the in- 
put state up( m \jri)u~ 1 will be mapped onto the output 
state u (g> up( out ^ (m)w _1 ® Another way of stating 
this condition is obtained by noting that if we express u 
as 



B. Proof of optimality 



Our proof of the optimality of the entangler (4.1) via 
the U-NOT gate is based on the recent idea of Gisin 
pl| , p2[ that the impossibility of instantaneous signal- 
ing generates upper bounds on the fidelity of particu- 
lar quantum-mechanical processes. To be more specific, 
we have shown earlier that the impossibility of the ideal 



u = exp(— i#e ■ B/2), 



(4.7) 



where e is a unit vector corresponding to the rotation 
axis and 9 is the rotation angle, then 



u(rh ■ B)u 1 = to' • <?, 



(4.8) 



where in' = R(e,9)rh. The rotation matrix, R(e,9), is 
the 3x3 matrix which rotates a vector about the axis e 
by an angle 9, and it is given explicitly by 







R(e,0) = exp(0e • K), 

where 


K, =|00-1 
1 

1 
K y =\ | , 
-10 

0-10 
A | 1 
v 

We have that 

up^ n \r7i)u- 1 = p {m) {Rm), 



(4.9) 



(4.10) 



(4.11) 



which will be mapped to p( out \Rm), so that the covari- 
ance condition can now be expressed as 

p {out) {Rm) = u®up (out) {m)u- 1 ®vT l . (4.12) 

Now let us examine the consequences of this relation. 
We shall first consider the terms linear in a and let R 
be a rotation about to by a very small angle 6. We have 
that 

a(Rfn) = Ra(m), (4.13) 
which for our choice of rotation becomes 

a(m) = (1 + 9m ■ K)a(m), (4-14) 



to • Ka(rh) = 0, 



(4.15) 



where to is a unit vector in the direction of to. This im- 
plies that to x a(m) — 0, so that a(rh) is parallel to to, 
and we can write a(m) — a(m) m. If we now substitute 
this result back into Eq. ( 4.13 ) and consider a general 
rotation i?, we have that 



a(Rm) = a(m). 



(4.16) 



This implies that o(to) is a constant, which, following 
pjj , we shall denote by rj. 

Now let us see what covariance implies about the terms 
quad ratic in a. Application of the covariance condition, 
Eq. ( 4.12 ), to these terms gives 



tj k (Rm) 



j',k' 



Rjj>Rkk'tj'k'(rn)- 



(4.17) 



If we again choose R to be a rotation about to by a small 
angle 9, we find the condition 



= ^)(m • K) jf t fk (m) + ^(m • K) hh rt jk ,{m). (4.18) 



k' 



If we choose to to be in the z direction, in particular 
to = we find, as did Gisin, that t xx 



and t xz — tzx tyz 
evaluated at to = z. 
signaling condition 



tyyi txy — ^yxi 

= t zy = 0, where all of these are 
We now want to impose the no 



p (out \z) + p (out \-z) = p^ out \x)+ p^ out \-x), (4.19) 



and to do so we need to find all of the density matrixes 
in the above equation in terms of tjk(z). This c an b e 
done by applying the covariance condition, Eq. (4.12), 
to p( out \z) and making the proper c hoice of R. When 



these results are substituted into Eq. ( 4.19 ) we find that 



txx(z) = tyy(z) — t zz (z), and we shall designate this 
common value by t(z). We then have that 



,(«•*)(!) = 



1 + 2r] + t 








1 - 1 
o 





2(t + it xy ) 
1 - t 



1 - 2?/ + t 
(4.20) 



The basis in which the matrix is expressed is {| + 
I, +z), | + z, —5), | — z, +z), | — z, —z)}, where a z \ ± I) = 
±|z). This matrix must be positive, which implies that 
the eigenvalues 

\(l±2v + t); i(l-*±2 v /*2 + t 2 w) (4.21) 

must be nonnegative. 

For an input state p^"^ (z) our desired output state is 
(| + z,—z) + | — z,+z))/-\/2, and this implies that the 
fidelity of p(° ut ) is 



T 



l + t 



(4.22) 



This is clearly maximized when t is as large as possi- 
ble, and examining the eigenvalues of p( out \ this happens 
when t xy = and t = 1/3. Substituting this into the ex- 
pression for the fidelity, we see that the maximum fidelity 
is 1/3. This means that the no-signaling constraint spec- 
ifies the upper bound on the fidelity of the symmetric 
entangling which is exactly the same one as achieved by 
the U-NOT gate. This proves that the entangling via the 
U-NOT gate is optimal. 



C. Remark 

We note that using the universal NOT gate one can 
also produce an entangled state of the form ( |3.l| ) . Specif- 
ically, the U-NOT gate allows Charlie (C) to produce an 
entangled state, consisting of l^) and one of two known 
states, which is shared by Alice (A) and Bob (B). In or- 
der to see how this can be accomplished it is usefu l to 
express the state on the right-hand side of Eq. (4.3) as 
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-(|*} A |$_) B c + WJ^-hc), 



(4.23) 



Science Foundation under grant PHY-9970507 and by the 
1ST project EQUIP under the contract IST-1999-11053. 



where 

|$ ) = = (|o)|i)-|i)|o)) 



V2 



V2 



is the singlet state. Charlie now measures his particle 
along the axis corresponding to the states |0) and |1). 
Whatever result he obtains for his particle, the other two 
particles will be in an entangled state shared by Alice 
and Bob. For example, if Charlie finds his partic le i n 
the state |1), Alice and Bob share the state in Eq. (3.1). 
Note that Charlie can choose the states with which the 
state |v?) will be entangled by choosing the axis along 
which to measure his particle. 

This implies that if we want to produce either the en- 
tangled state of |W) with |0) or the entangled state of \^>) 
with |1), and we don't care which one we get, this can 
be done with perfect fidelity. Perhaps a better way of 
stating this is that if we want to entangle \^f) with one 
of two orthogonal states, this can be done perfectly, and 
we will know with which state it is entangled. 



V. CONCLUSIONS 

In this paper we have studied various possibilities for 
entangling two qubits so the initial information about 
their preparation is preserved. We have studied a spe- 
cific situation when the state of one of the qubits is known 
while the second state is arbitrary. We have shown that 
entanglement via symmetrization in this case can be per- 
formed with a very high fidelity (much higher than the 
fidelity of estimation). This type of entanglement can 
be very useful for stabilization of the storage of an (un- 
known) quantum state of one qubit against environmen- 
tal interaction and a random imprecision p0| . We have 
shown that the U-NOT gate optimally implements the 
entanglement transformation |*) — * 1^)1*^) + 1*^)1*). 



This means that the transformation (4.3) is very special 
indeed - it describes the optimal cloning, the optimal U- 
NOT transformation as well as the optimal entangler. 



ACKNOWLEDGMENTS 

We thank Nicolas Gisin and Christoph Simon for help- 
ful discussions. This work was supported by the National 



[1] A. Peres: Quantum Theory: Concepts and Methods 
(Kluwer Academic Publishers, Dordrecht, 1993). 

[2] L. Hardy, Phys. Rev. Lett. 68, 2981 (1992). 

[3] L. Hardy Phys. Lett. A 167, 17 (1992). 

[4] C.H. Bennett, G. Brassard, C. Crepau, R. Jozsa, 
A. Peres, and W.K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 

[5] C.H. Bennett, and S.J. Wiesner, Phys. Rev. Lett. 69, 
2881 (1992). 

[6] A.M. Steane, Rept. Prog. Phys. 61, 117 (1998). 
[7] S. Massar and S. Popescu, Phys. Rev. Lett. 74, 1259 
(1995). 

[8] R. Derka, V. Buzek, and A. Ekert, Phys. Rev. Lett. 80, 
1571 (1998). 

[9] V. Buzek, M. Hillery, and R.F. Werner, Phys. Rev. A 60, 
R2626 (1999). 

[10] A. Barenco, A. Berthiaume, D. Deutsch, A. Ekert, 
R. Jozsa, and C. Macchiavello, SIAM Journal of Com- 
puting 26, 1541 (1997). 

[11] V. Vedral, M.M. Plenio, M.A. Ripin, and PL. Knight, 
Phys. Rev. Lett. 78, 2275 (1997). 

[12] N. Linden and S. Popescu, Phys. Rev. A 56, R3319 
(1997). 

[13] S. Hill and W.K. Wootters Phys. Rev. Lett. 78, 5022 

(1997) ; W.K. Wootters Phys. Rev. Lett. 80, 2245 (1998). 
[14] G. Alber, "Entangl ement and the line arity of quantum 

mechanics." e-print quant-ph/990110^. 
[15] D. Bures, Trans. Am. Math. Soc. 135, 199 (1969); 
A. Uhlmann, Rep. Math. Phys. 9, 273 (1976); 24 229 
(1986); see also W.K. Wootters, Phys. Rev. D 23, 357 
(1981). 

[16] A. Peres, Phys. Rev. Lett. 77, 1413 (1996). 

[17] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. 

Lett. A 223, 1 (1996). 
[18] N. Gisin and S. Popescu, Phys. Rev. Lett. 83, 432 (1999). 
[19] V. Buzek and M. Hillery, Phys. Rev. Lett. 81, 5003 

(1998) . 

[20] N. Gisin and S. Massar, Phys. Rev. Lett. 79, 2153 (1997). 
[21] N. Gisin, Phys. Lett. A 242, 1 (1998). 
[22] N. Gisin, private communication 



11 



